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We investigate ground-state properties and quantum phase transitions in the one-dimensional 
S = 1 spin-orbital model relevant to cubic vanadates. Using the density matrix renormalization 
group, we compute the ground-state energy, the magnetization and the correlation functions for 
different values of the Hund's coupling Jh and the external magnetic field. It is found that the 
magnetization jumps at a certain critical field, which is a hallmark of the field-induced first-order 
phase transition. The phase transition driven by Jh is also of first order. We also consider how the 
lattice-induced ferro-type interaction between orbitals modifies the phase diagram, and discuss the 
results in a context of the first-order phase transition observed in YVO3 at 77K. 

PACS numbers: 75.30.Et, 73.43.Nq, 75.10.Pq, 71.27.+a 



I. INTRODUCTION 

Orbital fluctuations in correlated electron systems have 
attracted much interest recently. They play a particu- 
larly important role in transition metal oxides, where the 
low-energy physics is controlled by the orbital as well as 
spin degrees of freedom. Orbitals of <i-electrons possess 
characteristic spatial anisotropy, which sometimes results 
in remarkable phenomena such as the reduction of effec- 
tive dimensions. For instance, in cubic vanadates such as 
YVO3 and LaV03, an effective one-dimensional (ID) sys- 
tem is formed along the c-axis at low temperatures, which 
has been shown by a number of experimentsi2iMiSiSil on 
spin and orbital orderings. For YVO3, there are two 
magnetic phases: C-type antiferromagnetic order (ferro- 
magnetic along c-axis and antiferromagnetic in ab plane) 
at higher temperatures 77 K < T < 114K, and G-type 
spin order (antiferromagnetic in all three directions) at 
lower temperatures T < 77K^ As for the orbital sector, 
the G-type orbital-antiferromagnetic correlations develop 
below T ~ 200-ftT and it changes to the G-type order 
(in-phase alignment of orbitals along the c-axis) below 
T = 77 K. On the other hand, in LaV03, the magnetic 
order is always of the G-type and the orbital correla- 
tions are of the G-type^ On empirical grounds, difference 
in the ground state of these two compounds is usually 
attributed to larger lattice distortions in YVOsjS. More 
specifically, an increase in octahedral tilting in YVO3 
supports ferro-type alignment of orbitals along the c-axis, 
hence stabilizing the G-type spin order in the ground 
state of YVO3. This argument has been supported by 
electronic structure calculations. 8 In vanadates, there are 
two electrons in the threefold degenerate ti g level, and it 
is believed that orbital degrees of freedom are crucial for 
understanding the interplay between different orderings 
in these compounds and their highly unusual magnetic 
propertiesi 2 : 3 : 4 : 5 : 6 : 7 

Based on the above experiments, an effective spin- 



orbital model has been proposed, 9 which may properly 
describe low-energy physics of the cubic vanadates. The 
ID version of this model consists of the S = 1 Heisenberg 
spin model with two- fold degenerate orbitals [see Eqs.QJ 
and It has been demonstrated that this effective ID 
model contains rich and essential properties inherent in 
spin-orbital systems^ A similar quasi- ID orbital model 
has been also proposed for LaV03. n 

In this paper, by calculating the ground-state energy, 
the spin and orbital correlation functions, and the magne- 
tization curve by means of the density matrix renormal- 
ization group (DMRG) method^ we investigate quan- 
tum phase transitions of the ID 5=1 spin-orbital model, 
which may be relevant for explaining the magnetic prop- 
erties of cubic vanadates, in particular, for understand- 
ing the phase transitions observed in YVO3. We ob- 
tain the phase diagram in the plane of the Hund's cou- 
pling Jh and the external magnetic field h. To make 
the model more realistic, we introduce also an additional 
ferromagnetic interaction V to the orbital sector, which 
may drive the system to an orbital-ferromagnetic phase. 
We discuss the importance of this interaction for the low- 
temperature phase of YVO3 below 77K. 

This paper is organized as follows. After brief dis- 
cussion of the model in Sec. II, we present in Sec. 
Ill the DMRG results for the ground-state energy, the 
spin/orbital correlation functions in a magnetic field, and 
show the magnetic phase diagram at absolute zero. In 
Sec. IV, we address the effects of the ferromagnetic or- 
bital coupling driven by lattice distortions. Brief sum- 
mary is given in Sec. V. 



II. MODEL 

We start with ID version of an effective Hamiltonian, 
which has been proposed to describe the spin and orbital 
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properties of cubic vanadates&iS 

"H = J y^[p(S t ■ Sj+i + 1) Jj,j+i + Kj,j+i] 



III. MAGNETIC PHASE DIAGRAM 



(1) 



where Si is the 5 = 1 spin operator at the i-th site. </(= 
4t 2 /U) is the superexchange interaction scale, which we 
will take as the energy unit in the following discussions. 
Jij and Ki j are the operators acting on the doubly- 
degenerate orbital degrees of freedom, and are given as 



jij = (1 + 2R) ( Ti ■ T 

1 

4 



K itj = R [Ti-Tj 



r [ T*T* 



-R, 



(2) 



where is the T = 1/2 pseudo-spin operator for the 
orbital sector. The terms proportional to R = 77/(1 — 3?y) 
and r = 77/(1 + 277) with 77 = Jh/U originate from 
the Hund's coupling (normalized by the on-site Coulomb 
interaction U). The effective model Q is derived via 
the second-order perturbation calculation in t/U (t is 
the electron hopping amplitude) for the cubic vanadates, 
where two correlated electrons occupy 7j 2g orbitals at each 
site. The microscopic derivation can be found in Ref. 9. 
Physically, the external field h in the above Hamiltonian 
can be regarded as the interchain molecular field present 
in a realistic 3D system^ 

For h = r\ = 0, it has been concluded that the 
ground state of the model JTJ is the orbital valence bond 
(OVB) solid state for which both of the spin and or- 
bital sectors are in a dimer-type singlet state with gap- 
ful excitationsAi2ii4 Furthermore, it was obtainec&i^ 
that a first-order transition from the OVB solid state to 
the fully polarized spin state with gapless orbital singlet 
may occur around a critical value ?7 C ~ 0.11. Since the 
analysis of Refs. ItWIoI is based on the finite-temperature 
DMRG method, it is difficult to discuss the behavior at 
lower temperatures. It is thus desirable to study zero- 
temperature properties to establish the precise phase dia- 
gram, which is addressed below in this paper. We further 
investigate the effects of a magnetic field, as well as of an 
additional ferromagnetic orbital coupling (which will be 
specified in Sec. IV). 

Before closing this section, we make a brief com- 
ment on the spin-orbital model. The 5 = 1 model 
may be regarded as a specific extension of the 
5 = 1/2 spin-orbital model, which has been studied 
extensively 15 i 16 i 17 i 18 ' 19 i 20 i 21 The ground-state phase di- 
agram of the latter model has been already established 
by the DMRG method^V ! 21 A simple isotropic 5=1 
generalization leads to the model which has the OVB 
ground statei 10 ' 14 We will see below that the extended 
model, which includes anisotropic orbital interactions, 
the magnetic field and also an additional ferromagnetic 
orbital interaction, exhibits interesting quantum phase 
transitions. 
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FIG. 1: The energy per site calculated by the infinite DMRG 
method as a function of the normalized Hund's coupling 77: 
the OVB solid state (triangle) and the spin-ferromagnetic 
state (circle). The inset shows the detail of the 77-dependence 
in the vicinity of r\ c . 
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FIG. 2: Spin (upper panel) and orbital (lower panel) corre- 
lation functions as a function of the Hund's coupling r\. The 
site index i is even (solid line) or odd (dashed line) . The even- 
odd i-dependence reflects a dimer property of the OVB solid 
state. 

We first consider the effect of the Hund's coupling. To 
this end, we precisely determine the critical value 77,, be- 
tween the OVB solid state and the spin-ferromagnetic 
state. We show the 77-dependence of the energy and the 
correlation functions calculated by the infinite DMRG 
method 12 for ft. = in Figs. □ and El It is seen from 
Fig. n that for small values of the Hund's coupling, the 
energy of the OVB solid state is lower than that of the 
spin-ferromagnetic state. In this region, both of the 
spin and orbital correlation functions show the even-odd 
i-dependence reflecting a dimer character of the OVB 
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solid state (see Fig. |2J). More precisely, for i = even, 
two adjacent orbitals form nearly perfect dimer pairs 
(namely, (T^ • TV|_i) is close to —0.75) while the spins 
prefer the ferromagnetic alignment with a positive value 
of (Si ■ Si+i) ~ 0.9. On the other hand, for i = odd, 
the value of the orbital correlation is almost zero while 
spins form a singlet pair with (Sj • Sj_|_i) ~ —1.5, which 
is much less than the value (= —2) expected for a spin-1 
isolated dimer £L The above behavior is characteristic of 
the OVB solid state. Increasing 77, the energy of the 
OVB solid state is getting closer to that of the spin- 
ferromagnetic state, and eventually a first-order phase 
transition occurs between these two phases at a critical 
value of r\ c = 0.1055 ± 0.0005. The detail around the 
transition is shown in the inset of Fig. ^ Beyond the 
critical value r\ c , the even-odd i-dependence of the cor- 
relation functions disappears and both of the spin and 
orbital correlations become spatially uniform with 

(S?SI +1 ) = 1, (T i .T 4+1 ) = -ln2+i (3) 

as should be expected^ Our conclusion that the 
ground state changes from the OVB solid to the spin- 
ferromagnetic state, and also the value of r/ c at which the 
transition occurs, is in good agreement with Ref. llOL 

The effective spin (orbital) exchange coupling con- 
stants, which are determined by local orbital (spin) cor- 
relations, change drastically at the phase transition. Just 
below rj c , the spin interactions J s — (Jj i j+i)/2 are equal 
to Jp ~ —0.38 within the ferromagnetic bonds, and 
Jaf — 0.05 between the spins belonging to different 
dimers. Above the transition, J s ~ —0.21 in all the 
bonds, supporting uniform Heisenberg-like orbital dy- 
namics. Similarly, the orbital exchange constant J or b = 
(1 + 2-R)(S; • Si+i + l)/2 + R changes from: J orb ~ 1.4 
(~ 0.02) in strong (weak) bonds below r/ c , to the uniform 
value J or b — 1.46 above r\ c . 

Let us next consider the magnetic properties of the 
model |J*EJ for finite values of h. We have calculated 
the magnetization curve as well as the correlation func- 
tions. The results of the magnetization are shown in 
Fig. El for two typical choices of the Hund's coupling. 
As seen from the case of 77 = in Fig. [3J the ground 
state in small fields is still in the OVB solid state, ren- 
dering the magnetization zero up to a critical field h c % 
due to the existence of the spin gap. Beyond the criti- 
cal field hcii the spin sector becomes gapless, resulting 
in the square-root increase of the magnetization near the 
critical field h c i- We note that this transition in the spin 
sector is a kind of Pokrovsky-Talapov transition, i.e. an 
insulator to gapless-liquid transition triggered by filling- 
control.^ 4 Therefore, the ferro/antiferromagnetic alter- 
nating spin correlation still exists even in the gapless 
phase. The magnetization then smoothly increases un- 
til we encounter a sudden jump at the second critical 
field h C 2, where the system undergoes a first-order tran- 
sition to the fully polarized spin state with the orbital 
sector forming a gapless T = 1/2 Heisenberg pseudo-spin 



chain. Note that in contrast to the spin sector, the or- 
bital sector is always gapful for < h < h C 2. For n = 
(77 = 0.05), the corresponding critical field is h C 2 = 0.122 
(h c2 = 0.065). 

Shown in Fig. 0] are the correlation functions for 
77 = under magnetic fields. As mentioned above, in 
weak fields, both of spin and orbital correlations exhibit 
the even-odd dependence characteristic of the OVB solid 
state. Their values are unchanged up to h c \, beyond 
which the even-odd dependence is gradually suppressed. 
Notice that the spin correlation function shows a cusp 
singularity at h c \ while the orbital correlation function 
smoothly changes. This is because the spin sector shows 
an insulator to gapless-liquid transition, whereas the or- 
bital sector still stays in the gapful phase, which is af- 
fected by magnetic fields indirectly via the spin-orbital 
interaction. In magnetic fields beyond h c i, the spin and 
orbital correlation functions take the values given in @ 
corresponding to the fully polarized spin state with a 
gapless orbital Heisenberg chain. Similar behavior in the 
field dependence of the correlation functions is observed 
for other choices of the Hund's coupling, as far as the 
ground state is in the OVB solid state at h = 0. 
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FIG. 3: Magnetization curves for r\ = (solid line) and 
77 = 0.05 (dashed line). The magnetization becomes finite 
for h > h c i and exhibits a jump at h C 2 characteristic of the 
first-order phase transition. 

We thus end up with the zero-temperature magnetic 
phase diagram, which is shown in Fig. |3J There are 
three distinct phases in the 77-/1 plane. In the region of 
weak fields and weak Hund's couplings, the OVB solid 
state is stabilized in the phases I and E; the orbital sec- 
tor is gapful while the spin sector is gapful (gapless) in 
the phase I (I) . The transition from the phase I to II is of 
Pokrovsky-Talapov type. On the other hand, the transi- 
tion from the phase II to the spin-ferromagnetic phase 11 
is of first order. 

For the spin/orbital dimerized state at 77 = 0, we es- 
timated the spin gap from the magnetization curve and 
obtained A s ~ 0.01J. This is actually much larger than 
the Haldane gap that one would expect from a simple 
picture of spin-2 dimers coupled by a weak antiferro- 
magnetic exchange. Indeed, at 77 = an effective spin 
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FIG. 4: Spin (upper panel) and orbital (lower panel) corre- 
lation functions for r\ = in a magnetic field. The site index 
i is even (solid line) or odd (dashed line). 



exchange constant between the dimers is estimated as 
J eff = l(T t -Tj + \) ~ 0.036. This would give the spin-2 
Haldane gap about 0.085 J eff ~ 0.003 only,S^ The large 
spin gap in the spin-orbital model is due to a dynamical 
coupling between two sectors, as discussed in Ref. 
We notice that our result for the spin gap is still much 
smaller than A s ~ 0.041 J, obtained in Ref. ITfl from the 
finite-temperature spin susceptibility \ computed by the 
transfer-matrix DMRG method. As the latter method is 
available at finite temperatures only, Ref. ^3 estimates 
the spin gap by extrapolating the data with a function 
of x ~ exp(-A s /T)/v / T. We think that the above dis- 
crepancy may arise from this (somewhat ambiguous) ex- 
trapolation procedure. 



IV. EFFECTS OF THE FERROMAGNETIC 
ORBITAL INTERACTION 

For representative values of Jjj ~ 0.68eV and U — 
^■Jh — 4.5eV for vanadates^ one obtains that the ra- 
tio 1] = Jr/U ~ 0.11 is slightly above its critical value. 
Therefore, our results are consistent with the experimen- 
tal fact that the ferromagnetic chain is stabilized along 
the c-axis in the ground state of LaV03- 6 However, the 
superexchange model alone is not sufficient to un- 
derstand the origin of the spin-G type ground state in 
YVO3, and also to understand the physics behind the 
spin/orbital reordering at T — 77 K. In order to describe 
low-temperature properties of YVO3 more precisely, it is 
desirable to take into account the pure-orbital interaction 
in addition to the spin/orbital couplings involved in the 
Hamiltonian (JTJ. The orbital interaction to be included 
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FIG. 5: Magnetic phase diagram: Phase I represents the 
OVB solid state with zero magnetization (spin-dimer state). 
Phase H is also in the OVB solid state, but with a finite mag- 
netization < (Sf) < 1. Phase EI is the spin- ferromagnetic 
phase with a gapless T = 1/2 orbital chain. The first-order 
phase transition occurs at the boundary denoted by the solid 
line. 

has the following Ising-likc fornifS, 

H = -VY,T?T?+v (4) 

i 

This describes a ferromagnetic orbital coupling, which 
competes with the antiferromagnetic orbital coupling in 
the Hamiltonian l[TJl. Physically, this term is important 
because apart from the superexchange mechanism there 
is also a lattice contribution to the orbital interactions. 
As mentioned in the introduction, tilting of octahedra 
in perovskites (caused by ionic mismatch effects) leads 
to a tendency that orbitals are aligned ferromagnetically 
along the c-axis. This effect is stronger when we go 
from La- to Y-based compound^ We thus expect that 
the ferromagnetic Ising orbital coupling along the c-axis 
may be particularly relevant for YVO3, in order to ex- 
plain why its ground state shows spin-antiferro/orbital- 
ferromagnetic structure along the c-axis. 

We choose here three different values of the Hund's 
coupling and calculate the energy as a function of V. 
In Fig. Eta), the ^-dependence of the energy for three 
competing states is shown for the OVB ground state 
(i] = 0.05). There is a direct first-order transition from 
the OVB state to the orbital-ferromagnetic state. On the 
other hand, when the Hund's coupling is slightly above 
the critical value r\ c (r\ = 0.11), we encounter different be- 
havior, as seen in Fig. Etb). Namely, the ground state is 
the spin-ferromagnetic (orbital-ferromagnetic) state for 
small V (large V), but for intermediate V, the OVB 
solid state is stabilized, which results from the compe- 
tition of the above two (spin-ferromagnetic or orbital 
ferromagnetic) interactions. Stabilization of the orbital 
disordered OVB state by finite V-interaction (induced by 
lattice distortion) is a remarkable result. The physics 
behind this observation is that the lattice-driven inter- 
action J3J introduces frustration into the orbital sec- 
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tor, competing with antiferro-type alignment of orbitals 
due to the superexchange process. For even larger 77 
[Fig. Etc)], V induces a direct phase transition from 
the spin-ferromagnetic state to the orbital-ferromagnetic 
state. 
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FIG. 6: The energy per site as a function of the orbital cou- 
pling V for several choices of 77; the OVB solid state (triangle), 
the spin-ferromagnetic state (circle), orbital-ferromagnetic 
state (square). Notice that the first-order transition occurs 
once for (a) and (c), while it occurs twice for (b) when the 
value of V is increased. 

The phase diagram thus obtained in the 77-V plane is 
shown in Fig. [7] There are three distinct phases in this 
figure. For small 77 and V, the OVB solid state is sta- 
bilized in the phase I, which is driven to the orbital- 
ferromagnetic phase IV with the increase of V. In 
the latter phase IV, the spin sector is described by an 
5=1 isotropic antiferromagnetic spin chain, and thus 
its ground state is a spin liquid with Haldane-gap excita- 
tions. In the spin-ferromagnetic phase IE, the orbital sec- 
tor is described by the XXZ pseudo-spin model, which 
results from (JTJ and Q: 



H = J^,T i -T i+1 -Vj2 T l T i 



i+i- 



(5) 



Here J = J/(l — 377), and we have discarded irrelevant 
constant terms. According to the well known results for 
the 5=1/2 XXZ spin chain with ferromagnetic Ising 
anisotropy, 2 ^ the ground state is in a gapless liquid phase 
for V < 2 J, and in a fully polarized ferromagnetic phase 
for V > 2J. Therefore it might be possible to have the 
orbital- ferro/spin- ferromagnetic ground state in the re- 
gion V > 2/(1 — 377) in units of J. However, it turns out 
that since the first-order transition from the phase IH to 
IV occurs at the phase boundary shown in Fig. [7| there 
is no available window of the parameters which can sta- 
bilize the orbital- ferro/spin- ferromagnetic ground state. 

The ferromagnetic orbital interaction V may be rel- 
evant to the first-order phase transition observed in 
YVO3 at T = 77 K, where the spin ordering changes 
from the G-type (below 77K) to the C-type (above 
77A") structure. If we focus on the pattern of order- 
ing in the c-direction, the transition occurs from the 



spin-antiferro/orbital- ferromagnetic ground state to the 
spin-fcrro / orbital-antifcrromagnetic state. Furthermore, 
recent neutron scattering experiments have revealed a 
substantial modulation of the ferromagnetic spin cou- 
plings along the c-direction at T > 77 K, which indicates 
the appearance of orbital dimer correlations at finite 
temperatures^ According to the results obtained here, 
such a temperature-driven first-order transition may pos- 
sibly occur when the system is located in the phase IV 
near the tricritical point in the phase diagram (Fig. EJ, 
and thus close to the phases IH and I. For example, 
V ~ 1.5 may be a reasonable choice for 77 ~ 0.12, which 
may give rise to the competition of the phases at finite 
temperatures. This is in fact quite realistic situation, as 
the atomic value of Hund's coupling parameter for vana- 
dates is indeed close to this value Finally, we wish 
to mention that the gapless orbital-liquid (gapful spin- 
liquid) state in the phase II (IV) is expected to show 
an orbital (spin) antiferromagnetic order when realistic 
three-dimensional effects are properly taken into account, 
being consistent with experimental findings. 
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FIG. 7: Phase diagram in the 77-V plane; Phase I is the OVB 
phase, II is the spin-ferromagnetic phase, and IV is the orbital- 
ferromagnetic phase stabilized by the orbital ferromagnetic 
Ising term. All the phase transitions are of first-order. Note 
that the orbital (spin) sector in the phase HI (IV) is a gapless- 
orbital (gapful-spin) liquid. 



V. SUMMARY 

We have investigated quantum phase transitions of the 
ID 5 = 1 spin-orbital model relevant to cubic vana- 
dates. By means of the DMRG calculation, we have de- 
termined the zero-temperature phase diagram in a mag- 
netic field. We have found that the magnetic field induces 
the phase transitions twice for the case of OVB ground 
state. The lower-field transition is of Pokrovsky-Talapov 
type, which drives the spin sector to a gapless spin liquid 
while keeping the orbital sector always gapful. On the 
other hand, the high-field transition is first order to the 
spin- ferromagnetic state, where the orbital state is driven 
to a gapless isotropic T — 1/2 orbital chain. 
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We have considered the evolution of spin and orbital 
correlations as a function of the Hund's coupling, and ob- 
tained that the ground state of the superexchange model 
is spin- ferromagnetic for the realistic values of Jr- This 
is consistent with the formation of ferromagnetic spin 
chains along the c-axis in LaVC>3. 

The introduction of the ferromagnetic Ising coupling 
V to the orbital sector induces a first-order phase tran- 
sition to the orbital- ferro/spin-antiferromagnetic state. 
This may explain the stabilization of the G-type spin or- 
der observed in YVO3 at temperatures below 77K. We 
found also that the intermediate values of U-interaction 
may stabilize the spin/orbital dimer phase. This phe- 
nomenon reflects competition between the superexchange 
and lattice effects, and, surprisingly enough, it occurs for 
Jff- values realistic for vanadates. The proximity of differ- 
ent spin/orbital states in the ground state may result in 
an entropy-driven phase transition at finite temperature, 



similar to that observed in YVO3. In the future study, it 
is desirable to discuss whether such a tcmpcraturc-driven 
first-order transition may be indeed possible in the spin- 
orbital model including the ferromagnetic orbital inter- 
action and also the interchain coupling effects. This may 
provide a key to clarify some essential spin/orbital prop- 
erties of YVO3 at low temperatures. 
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